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NOTATION  (Cont’d) 


=  dw/dr  j  radial  slope  of  a  deflected  element  of  the  plate 
initial  radial  slope  of  the  plate 
function  of  p  defined  by  Eq. (2.2.10) 


INTROEUCTION 


Heat  energy  may  cause  the  failure  of  a  structure  in  a  number  of 
ways;  in  particular,  thermal  stresses  may  cause  instability  of  the 
configuration  of  a  structural  element. 

This  report  summarizes  the  methods  and  conclusions  of  an  investi¬ 
gation  on  thermally  induced  instability  of  annular  plates.  The  research 
work  here  described  is  limited  to  cases  in  which  the  temperature  veurl- 
atlon  is  symmetric  about  the  axis  of  the  annular  plate.  Moreover,  the 
temperature  is  assumed  to  be  uniform  across  the  thickness  of  the  plate. 
Additional  siiiQ)llfying  assumptions  will  be  listed  later. 

In  the  first  section  of  this  report,  consideration  is  given  to 
the  problem  of  determining  the  radish,  temperature  distribution  corre¬ 
sponding  to  prescribed  boundary  conditions. 

In  the  second  section  the  problem  of  determining  theoretically  the 
occurrence  of  instability  for  a  given  plate  and  a  given  transient  temper¬ 
ature  distribution  is  analyzed  and  a  method  is  presented  to  obtain 
approximate  solutions.  The  method  is  applied  to  two  different  cases: 
in  the  first  the  temperature  distribution  is  obtained  theoretically  and 
corresponds  to  idealized  heat  transfer  conditions;  and  in  the  second 
the  ten^erature  distribution  is  assumed  to  be  known  from  em  experimental 
Investl^tion. 

The  effect  of  initial  imperfections  in  the  flatness  of  the  plate 
is  considered  and  a  Southwell  plot  is  used  to  determine  the  critlced. 
temperature  from  experimental  data  when  the  initial  imperfections  are 
smeill  and  axlsymmetrlc . 

The  third  section  is  a  summary  of  the  experimented  methods  applied 
to  the  collected  data  on  temperature  euid  deflections.  A  description  of 
the  tests  performed  is  presented  euid  the  results  eure  compeured  with  the 
theoretical  predictions. 

Finally,  eui  appendix  covers  the  interpretation  of  the  interference 
fringes  that  are  obtained  by  the  Moire  technique. 
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SECTION  I 


THE  AXISYMMETRIC  PROBLEM  OF  HEAT  TRANSFER 
IN  A  CIRCULAR  PLATE 


1.1  General  Problem 

In  the  general  case  the  plate  under  consideration  nay  receive  beat 
energy  through  conduction  from  an  internal  heat  source  or  through  con¬ 
vection  or  radiation  from  an  external  source.  In  either  case  it  is 
assiimed  that  the  heat  input  through  the  boundary  of  the  plate  is  axl- 
symmetric.  During  the  transient  period  the  heat  flows  throu^  the  plate 
in  accordance  with  the  laws  of  conduction  and  is  partially  stored  In 
the  plate  and  partially  dissipated  to  the  environment  by  convection  and 
radiation.  When  the  steady  state  is  reached  the  heat  dissipation  to 
the  environment  balances  the  heat  input  emd  the  heat  stored  remains 
constant. 

The  general  problem,  as  stated  above,  is  nonlinear  and  difficult 
to  solve.  No  attempt  will  be  made  here  to  solve  it.  Instead,  slnqpll- 
fying  assumptions  will  be  Introduced  in  order  to  linearize  the  problem. 

1.2  Heat  Transfer  without  Dissipation  to  the  Environment 

When  the  assumption  is  made  that  no  heat  transfer  occurs  from  the 
lateral  surfaces  of  the  plate  to  the  environment,  the  problem  is  linear 
and  is  analogous  to  that  of  em  infinite  cylinder. 

Solutions  to  this  problem  for  different  boundary  conditions  are 
readily  available  (see  Ref.l).  A  particular  case  of  interest  will  be 
treated  here. 

1.2.1  Annular  Plate  with  Constemt  Flux  at  the  Inner  Edge 
emd  Zero  Temperature  at  the  Outer  Edge 

The  solution  to  the  case  of  constant  flux,  F  ,  at  tbe 
inner  edge  (  r  «  a  ),  euid  zero  teogperature  at  the  outer  edge  (  r  ■  b  ) 
is:  . 
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vhere  the  are  the  positive  roots  of  the  equation 


Jj^(aa)y^(ba)  -  Y^(aa)j^(ha)  -  0 


(1.2.1.2) 


.  In  the  experimental  part  of  this  reseaireh  a  nuniber  of  tests  ' 
vere  conducted  for  an  annular  plate  specimen  of  the  foUoving 
characterlstios: 

Material:  mild  steel 

Inner  radius:  a  =  1.5  In* 

(1.2.1.3) 

Outer  radius:  b  ■  5  In. 

Dlffuslvlty:  k  •  0.12  cnP/tee 

With  these  particular  values  the  first  three  roots  of  equation 
(1.2. 1.2)  are; 

■  0.2194  cm”^ 

=  0.5507  cm"^  (1.2.1.4) 

0^  ■  0.8966  ca"^ 

After  substitution  of  these  values  in  (1.2. l.l)  the  expression  (t/F^)K 
can  be  plotted  versus  r  for  particular  values  of  the  total  heating 
tine  t  .  In  Flg.l  this  has  been  done  for  t  ■  50^  100,  200,  400  sec 
and  for  t  ■  (steady  state). 

For  values  of  tine  lower  than  50  seconds,  more  than  three 
roots  would  be  required  In  the  calculations.  If  0  • 
p  s  r/b  are  Introduced,  where  Is  the  temperature  at  the  Inner 
edge  (  r  ■  a  )  a  nondlmenslonal  plot  of  the  temperature  distribution 
for  the  peurtlculeu:  plate  under  consideration  is  obtained  (Fig. 2).  It 
is  interesting. to  note  that  9  is  Independent  of  the  heat  flvuc  F^  . 
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1.3  Heat  Tremsfer  vlth  Dissipation  to.  the  Environment. 

Experimental  Survey 

In  the  absence  of  a  theoretical  solution  accounting  for  the  dissi¬ 
pation  of  heat  to  the  environment,  an  experiment  was  performed  to 
estimate  the  order  of  magnitude  of  the  effect  on  the  temperature  dis¬ 
tribution.  An  annular  plate  of  the  dimensions  specified  In  (1.2. 1.3) 
vas  rigidly  supported  at  Its  outer  edge  by  means  of  a  massive  ring,  so 
as  to  reproduce  closely  the  botindary  condition  T»0  at  rab  of 
1.2.1.  An  electrical  heater  located  In  the  center  of  the  annular  plate 
provided  a  heat  source  of  reasonably  constant  heat  flux  at  r  »  a  . 

Figure  S  is  a  plot  of  the  tenqperature  variation  with  tise  for  a 
nundOer  of  points  on  the  plate.  The  Irregular  behavior  noticeable  between 
0  ai4  30  seconds  corresponds  to  the  period  which  the  heater  requires 
to  attain  steady  temperature;  .. 

Figure  is  a  plot  of  the  tenqperatUre  distribution  In  the  plate 
for  different  values  of  tlB».  This  is  again  plotted  nondlmenslonally 
in  Fig.5. 

The  tenqperature  distributions  for  t  >  200  from  Figs.  2  and  5 
are  superimposed  in  Fig.6  for  the  purpose  of  comparison.  The  coiqparison 
Is  particularly  significant  due  to  the  fact  that  the  nondimens lonal 
temperature  distribution  in  the  idealised  case  is  Independent  of  the 
heat  flux. 
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SECTION  II 


LINEARIZED  THEORY  OF  THE  AXISYMMETRICAL 
THERMAL  BUCKLING  OF  ANNULAR  PLATES 

2.1  Basic  Assumptions 

Small  deflections  are  as sinned,  as  Is  the  case  at  the  onset  of 
buckling,  and  consequently  the  usual  assumptions  of  the  theory  of  thin 
plates  with  small  deflections  will  be  made. 

Tenqierature  distribution  and  edge  conditions  are  perfectly  axl- 
symmetrlc.  The  annular  plate  is  perfectly  flat,  of  uniform  thickness 
at  a  uniform  temperature  and  free  from  stresses  In  its  initial  state. 
Temperature  variations  through  the  thickness  are  neglected.  The  material 
Is  perfectly  elastic,  homogeneous  and  isotropic. 

2.2  Differential  Equation  for  Axlsymmetrlcal  Bending 

The  equation  of  equilibrium  in  the  lateral  direction  of  an  element 
of  the  plate  in  cylindrical  coordinates  can  be  written  (Ref. 2) 

^  [j  ^  W]  -  -  §  (2.2.1) 

where  Q  denotes  the  shearing  force  per  unit  length  of  the  cylindrical 
section  of  radius  r  . 

With  the  assumption  of  small  deflections  and  In  the  absence  of 
lateral  loads,  the  shearing  force  satisfies  the  following  relation: 

Q=  {o^h)<p  (2.2.2) 

Substitution  in  (2.2.1)  yields 

which  Is  a  linear  homogeneous  second-order  equation  In  ^  is 

given  by  the  equation  (Ref. 3) 
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(2.2.4) 


a  =  -  aE 
r 


7/  7] 


where  c^,Cg  are  determined  by  the  bouiadary  conditions.  In  the  case 
of  an  annular  plate  with  fixed  outer  boundary  and  free  inner  boundary 
the  boundary  conditions  are: 


r=a 


(u)  ^  =  0 

'  'r»h 


(2.2.5) 


where  u  is  the  radial  displacement  given  by 


u  .  (1«)  O  i/  Trdr  +  Cj^: 


r  +  <=2  7 


(2.2.6) 


Substitution  yields:  ^  •'.•'V- 


a  =  -  2|  r.  r  Trdr  +  f  Trdrl 

^  V  Lb^(l-v)  ♦  a^(l+v)  "^a  ,  ^a  J 


If  the  nondimenslonal  parameters 


(2.2.7) 


p  =  r/b 

•  0(p)  -  T(p)/T„ 


(a/b  s  P  S  1) 


(2.2.8) 


are  Introduced,  vbere  is  the  ten^rature  at  r  ■  a  ,  and  are  sub¬ 
stituted  In  equations  (2.2.7)  9x0.  (2.2.5),  the  differential  equation 
(2.2.3)  can  be  written  in  the  fom 


fe[p  ^  ^ 


(2.2.9) 


-  6 


where 


t  = 


.(pV, 

I  (1-v)  +  a‘^(l+v)  "a/t  a/b 


(2.2.10) 


and 


-  12(l-v^) 

h 


(2.2.11) 


with  the  boundary  conditions 

■  0  (zero  slope  at  p  =  1  ) 

(2.2.12) 

(^  +  V  ■  0  .  (zero  moment  at  p  =  a/b  ) 

^'p»a/b 

anil 

^y^p:»l.*  I®  fe{p  Ip  ^  ®  (shear  force  at  p  -  1  ) 

This  last'boundary  condition  Is  identical  with  equation  (2.2.9)  for 

p  a  1  when  the  condition  (tp)  ,  =  0  is  applied.  This  is  due  to  the 

P** 

fact  that  equation  (2.2.9)  is  the  equation  of  equilibrium  In  the  lateral  . 

direction  of  an  element  of  the  plate.  The  condition  of  zero  shear  force 

at  the  inner  edge  is  also  satisfied  by  equation  (2.2.9)  &8  &  consequence 

of  condition  (2. 2. 9),  (o  )  ■  0  .  This  condition,  applied  to  (2.2.3) 

.  **  r»a 

results  in  (v)  >  0  . 

The  determination  of  ip  reduces  to  the  boundazy-value  problem  of 
solving  equation  (2.2.9)  vlth  the  boundary  conditions  (2.2.12). 

The  term  t  (2.2.10)  is  a  function  of  p  detennlned  by  the 
geometrical  and  elastic  characteristic  of  the  plate  and  the  te]iq>erature 
distribution  in  the  plate. 
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2.3  Determination  of  the  Critical  Temperatiure 

In  the  case  of  differential  equation  (2.2.9) ^  ^  =  0  is  always  a 
solution  of  the  problem.  Solutions  with  (p  4  ^  can  be  found  only  for 
discrete  values  of  the  parameter  .  Physically  this  means  that  the 
plate  may  adopt  different  modes  of  deflection  for  corresponding  differ¬ 
ent  discrete  values  of  the  tecqperature  .  The  least  of  these  values 

will  be  taken  as  the  critical  temperature-  T  ...  It  can  be  shown 

crit 

from  energy  consideration  that 


If  Tq  <  £  0  Is  the  stable  state  of  equilibrium 

If  >  <;?  E  0  is  unstable  and  some  solution  (p  ^  0  Is  stable. 

.  '  Exact  solutions  of  eqiiatlon  (2.2.9)  cannot  be  found  readily  for 
prescribed  temperature  distributions,  but  several  approximate  methods 
are  available  to  solve  ;the  problem.  ' 

■2.U  Application  of  the  Stodola-Vtanello  Method  to  the 
Determination  of  the  Critical  Temperature 

The  Stodola*-Vlanello  method  (see  Ret»U)  Is  a  method  of  successive 
approximations  which  gives  a  sequence  of  functions  that  converge 
precisely  to  the  first  , characteristic  value. 

Two  successive.' Integrations  of  equation  (2.2.9)  yield  the  expression 


(2.4.1) 


The  method  consists  of  starting  with  a  conveniently  chosen  first 
approximation  (p^  to  <p  In  the  right  side  of  (2.4.1).  After  the  two 
successive  Integrations  are  completed,  c^  and  c^  are  determined  in 
such  a  manner  that  the  expression  within  the  brackets  satisfies  the 
boundary  conditions.  This  will  result  in 


^2  “  Vl^P^ 


(2.4.2) 
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If  the  approximation  <p^  had  been  the  characteristic  function 
then  ^2  =  <p^  and  the  corresponding  characteristic  value  of  t  would 
be  given  by 


T 


O 


^l(p) 


(2.U.3) 


The  first  approximation,  will  in  general  not  give  a  con¬ 

stant  ratio.  If  ^2  is  taken  as  the  improved  approximation  and  the 
cycle  repeated,  it  can  be  shown  that  the  ratio  (p^{p)/f^{p)  tends  to 
a  constant  value  over  the  interval  (p^^l)  ^8  is  increased;  and 

that  this  constant  value  is  the  lowest  characteristic  value  t  . 


Successive  estimates  of  the  characteristic  value  may  be  obtained 

after  each  cycle  by  requiring  that  the  functions  p  and  <p  ,  =  t  f 

n  ’^n+l  o  n 

agree  as  well  as  possible  (in  some  sense)  over  the  interval.  One 
possible  formulation  is 


{2.k.k) 


The  method  cein  be  carried  out  either  algebraically  or  graphically. 
However,  difficulties  are  encountered  in  our  particular  case  when  an 
algebraical  treatment  is  attempted.  The  temperature  distribution  d 
is,  in  most  cases,  not  easily  expressible  mathematically,  and  even  when 
this  is  done,  it  usually  results  in  complicated  expressions  for  f 
which  make  the  integrations  cumbersome.  It  is  preferable  therefore  to 
carry  out  the  work  graphically  (see  Fig.9)« 

From  a  plot  of  the  nondimenslonal  ten^rature  distribution  6 
versus  the  nondimenslonal  radius  p  ,  function  )|r  is  computed  numerically. 
Next,  the  first  approximation  <p^  to  the  slope  function  is  selected. 

The  rapidity  of  the  convergence  of  the  method  depends  on  this  selection. 
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The  simplest  algebraic  expression  that  satisfies  the  boundary 
conditions  (2.2.12)  is 


2  fv/p„  -  P„(2+v)-|  ,  ^ 

’’I'f  *L  (i«)  -  v/p„  J 


(a.ii.s) 


This  function  gives  a  very  fast  convergence  since  It  is  a  good 
approximation  to  the  actual  slope  function  (see  Section  III,  Fig. 33). 

Next  the  product  iv*},  computed  numerically  and  plotted.  The 
first  Integration  is  performed  and  the  product 


r” 

pTl  «  P  J 


(2.U.6) 


is  plotted.  In  the  following  step  the  second  Integration  is  ceurrled  out, 
the  result  divided  by  p  ,  and  the  boundaiy  conditions  are  applied  to 
determine  c^  and  Cg  .  -  This  yields  .  ^ 


(2.4.7) 


Next  f^(p)  is  plotted  and  the  first  approximation  to  the  characteristic 


value  T  is  obtained  fr(» 

o  •  •  •  . 


■  .  I 

.  .  _ 

-  /... 


(2.4.8) 


and  the  characteristic  function  from 


^2  - 


(2.4.9) 
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Because  the  first  approximation  Is  good,,  due  to  the  already 
mentioned  suitability  of  the  approximation  given  in  to  the 

problem  under  consideration,  one  cycle  of  the  calculations  suffices 
to  determine  an  acceptable  value  of  . 

The  critical  temperature  is  then  given  by  (2.2.11): 


(t  .)  =  T 

o  crlt  o 


12(l-v^)b^a 


(2.4.10) 


It  is  of  interest  to  remark  that  is  determined  only  by  the 
parameters  a,  b  and  the  Poisson  ratio,  plus  the  temperature  distri¬ 
bution  and  the  boundary  conditions,  that  is 

=  T^(a,b,v,6(p)  ,  boundary  conditions 


but  is  independent  of  the  thickness  of  the  plate.-  Equation  (2.4.10)  can 
be  written,  if  so  desired,  as 


1 

12(l.v=)b^a 


(2.4.11) 


2.4.1  Remarks  on  the  Applicability  of  the  Method 

In  the  preceding  paragraph  a  method  was  discussed  that^ 
starting  from  a  particular  nondlmenslonal  temperature  distribution 
6(p)  ,  yields  the  corresponding  value  of  the  critical  edge  temperature. 

Most  practical  problems  of  beating  result  In  transient 
conditions  In  which  9  Is  a  function  of  time  6  ■  d(p,t)  .  If  this 
is  the  case,  successive  approximations  a:ra  necessary  to  determine  the 
critical  temperature.  The  method  given  In  2.4  can  be  applied  to  deter¬ 
mine  critical  temperature  values  for  a  nuinber  of  values  of  t  ,  so  that 
a  plot  of  T^^  versus  time  can  be  obtained.  If  the  edge  ten^erature 
T^(t)  versus  time  Is  plotted  on  the  same  graph  the  Intersection  of 

the  two  curves  yields  both  T  and  t  . 

cr  cr 


-  11 


2.5  Application  to  a  Particular  Idealized  Case:  Thermal  Buckling 
of  an  Annular  Plate  with  Constant  Heat  Flux  at  the  Inner  Free 
Edge,  Zero  Temperature  at  the  Fixed  Outer  Edge  and  no  Dissi¬ 
pation  of  Heat  to  the  Environment 

In  paragraph  1.2.1  the  temperature  distribution  corresponding  to  a 
particular  annular  plate,  characterized  by  the  pcurainetrlc  values  given 
In  (1.2. 1.3),  vas  obtained  and  the  results  plotted  In  Flgs.l  suid  2. 

It  was  found  that  the  temperature  is  proportional  to  the  heat  flux.  The 
following  relation  can  be  written 

T(p)=  f(p,t)F^  (2.5.1) 


and  In  particular,  for  the  temperature  at  the  Inner  edge 

T^  =  f(t)FQ  (Fig.8) 


(2.5.2) 


On  the  other  hand,  6  =  T/T^  (Fig. 2)  is  independent  of  F^  .  If 
we  assume  that  buckling  occurs  at  a  particular  time,  for  Instance 
t  =  50  sec,  we  cein  calculate  from  Fig. 2,  which  provides  6^  ,  the 
corresponding  critical  temperature  .  (T  ) 


0  crit 

From  Fig. 8  we  know  that  for  t  =  50  sec 


50 


(1000  t/T  )  =  2;686 

o't=50 

If  T  =  (T  )  we  obtain  the  value  of  F  for  which  t  =50  sec 

o  °  °  cr 

F  *  1000  (T  ) 

0  2.^86  '  0  gpit 


50 

This  procedure  has  been  used  when 

t  =  50,  100,  200,  400  sec,  and  t  =  * 


and  the  results  have  been  plotted  In  Fig.  9* 

The  thickness  of  the  plate  has  been  left  as  a  design  paraiaeter  In 
accordance  with  (2.4.11). 
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Some  interesting  conclusions  can  be  drawn  from  Fig. 9* 

(a)  There  is  a  particular  value  of  F^/h  ,  that  we  shall  call 

(f  /h^)  ,  such  that 

O'  crlt 

for  F7h  i  (F  /h  )  buckling  does  not  occur: 

o  0  crit 

/  2  '  2 

for  F  /h  >  (F  /h  )  buckling  occurs  after  a 

°  o  crit 

finite  time  t  =  t 

cr 

The  critical  value  (F  /h  )  is  obtained  by  assuming  that 

•  .  o  crit 

buckling  occurs  at  t  ■  •  and  by  taking  the  nondimensional 
temperature  distribution  in  the  steady  state  as  the  nondimen- 
slpnal  temperature  distribution  at  buckling. 

2 

The  corresponding  value  of  (T  /h  )  is  a  mininrum  and  can 

°  crit 

be  taken  as  a  lower  bound  on  the  heat  flux. 

(b)  For  F  /h^  >  (F./h^)  the  larger  F  /h^  ,  the  larger 

,  2  °  °  crit  °  ,  p 

(T./h  )  and  the  smaller  t  .  When  F  /h  tends  to  »  . 

o  crit  2  cr  o' 

t  -*  0  and  (T  /h  )  tends  to  »  . 

'  p'  'crit 

This  can  be  proved  by  considering  the  process  of  graphical 
calculation  which  leads  to  the  determination  of  (T  /h  ) 

crit 

If  it  is  assumed  that  buckling  occurs  after  an  infinitesimal 
period  of  tin«,  the  corresponding  nondimensional  temperature 
distribution  0  will  be  a  line  very  close  to  the  axis  of 
abscissas  and  the  area  under  this  line  will  be  infinitesimal. 


The  successive  integrations  lead  to  an  expression  of  the 
following  type: 


.  ^  ^  '^iP  +  *=2  =  *  Vl^P^ 

where  €(p)  is  an  infinitesimal.  If  the  boundary  conditions 

fj^(l)  -  0  Cj^  +  Cg  -  e 

Sfl  f 

■^+v~-0  Cg-Ac^ 
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are  applied,  one  obtains 


and 


f^Cp)  =  € 


Hence 


T  = 

o 


/  <P-,^ 


finite 


€dp 


0  . 


crlt 


Now,  from  Fig. 8  we  see  that 


11m  =  0 

t  ^  0  o 


Consequently 


11m  F  ■*  <*> 
t  -  0  ° 


p 

Figures  8  and  9  are  combined  In  Fig. 10,  where  (T  /h  )  and 

°  crlt 

p  p 

T^/h  ,  are  plotted  against  time  for  different  values  of  yJIel.  For 

any  given  value  of  F  /h^  ,  the  Intersection  of  the  corresponding  curve 

2  ^2 

T  /h  with  the  curve  (T  /h  )  determines  the  critical  time  and 

o  o  crlt 

the  critical  tenqjerature.  The  conclusions  drawn  under  (a)  emd  (b)  are 

2  2 
visible  from  Fig. 10.  For  values  of  F  /h  smaller  than  (F  /h  ) 

o'  'o'  'crlt 

there  Is  no  Intersection,  which  means  that  buckling  does  not  occur. 
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2.6  Application  to  Experiment:  Thermal  Buckling  of  ein  Annular 
Plate  Centrally  Heated  at  the  Inner  Free  Edge  and  Fixed 
at  the  Outer  Edge 

In  paragraph  1.3  a  summary  was  given  of  an  experimental  study  of 
the  temperatxire  distribution  over  a  p«u*tlcular  specimen  plate  under 
specified  heating  and  boxmdary  conditions.  The  details  of  the  experl 
meht  are  described  In  Section  III. 


Figures  3^  1;  and  3  summarize  the  experimental  data  obtained  regard¬ 
ing  the  ten^ierature  distribution  during  periods  of  time  ranging  from  0 
to  300  seconds;  zero  corresponds  to  the  time  when  the  heat  Is  applied. 

The  time  of  buckling  can  be  determined  in  a  first  approximation  by 
applying  the  method  described  In  2.k  and  2.h.l, 

Starting  from  the  values  of  9  plotted  in  Fig.5>  the  corresponding 

values  of  (T  /h  )  are  calculated  and  plotted  against  time  (Fig.ll). 

°  crlt 

Next  the  values  of  T  /h  measured  during  the  process  of  heating  are 

plotted  against  time  on  the  same  graph.  The  intersection  of  (T  /h  ) 

2  .  °  crlt 

and  T^/h  determines  the  critical  point 

In  this  part iculsur  case  the  critical  values  are 


t  =  148  sec 
cr 


The  thickness  of  the  plate  used  In  the  experiment  was 


h  =  0.058  In. 

With  this  value  one  obtains 

(T„)  =  119. 5°P 

”  crlt 

For  the  inirpose  of  comparison  the  curve  (T  /h  )  corresponding 

°  crlt 

to  the  Idealized  case  in  which  dissipation  to  the  environment  was  neg¬ 
lected^  Is  plotted  as  a  dashed  line  In  the  same  graph. 
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In  the 


In  article  2.6  it  waa  noted  that  the  curve  (T  /h^^ 

0  crit 

idealized  case  was  an  invariant  with  respect  to  the  heat  flux  at  the 
heat  source.  Whether  this  applies  to  practical  cases  with  dissipation 
to  the  environment,  is  a  problem  that  requires  further  investigation. 
Such  an  investigation  could  consist  of  an  attempt  to  determine  theoret¬ 
ically  the  tengierature  distribution  when  dissipation  to  the  environment 
is  taken  into  consideration,  or  an  extensive  series  of  tests  from  which 
data  on  the  temperature  distribution  for  different  values  of  the  heat 
fl\uc  could  be  obtained. 

In  article  2.6  it  was  found  that  a  critical  value  of  the  heat  flux 

F  and  a  lower  bovind  for  (T  /h  )  could  be  determined  from  con- 

°  °  crit 

slderation  of  the  steady-state  ten^jerature  distribution.  If  the  steady- 

state  temperature  distribution  for  the  case  considered  in  this  article 

2 

were  known,  an  asymptote  could  be  determined  for  the  curve  (T  /h  ) 

o  crit 

and  the  corresponding  critical  value  of  F  and  lower  bound  of 

(T  /h  )  could  be  found, 

o  crit 

2.7  Effect  of  Initial  Imperfections 

Let  us  assume  that  the  initial  deflection  of  the  plate  is  axl- 
symmetrlc,  is  small  compared  with  the  thickness  of  the  plate,  and  is 
denoted  by  w^(r)  ,  and  that  the  additional  deflection  is  denoted  by 
w(r)  . 

In  article  2.2,  equation  (2.2.1),  the  equation  of  equilibrium  of 
an  element  of  the  plate,  when  subjected  to  axlsynmietrlc  bending,  was 
given  in  cylindrical  coordinates  as 

dr  [r  dr  D 

The  left-side  member  of  this  equation  is  obtained  from  the  expres¬ 
sions  for  bending  moments  and  since  these  moments  do  not  depend  on  the 
total  curvature  but  only  on  the  change  in  curvature  of  the  plate,  the 
additional  deflection  w(r)  .  should  be  used  in  this  mesiber. 
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On  the  other  hand,  Q.  being  the  shear  force  per  unit  length 
normal  to  the  Ideal  flat  initial  plate,  and  not  to  the  actual  initial 
plate,  it  can  be  expressed  as 

Q  =  (p) 

Substitution  yields: 

^  (r<p)  =  -  ^  (v  ^2.7.1) 

This  is  a  non-homogeneous  differential  equation  which  has  solutions 

for  arbitrary  values  of  o  .  The  physical  reason  for  this  fact  is  that 

r 

deflections  will  appear  as  soon  as  stresses  develop  in  the  plate. 

If  it  is  assumed  that  is  proportional  to  ^  ,  Eq. (2.7.1)  can 
be  written 

This  leads  to 

T  ..  =  T  (1+c)  r,  T  =  T  (l  +  ~ ] 

crit  o'  '  o\  <p  J  o^  w^y 

or 


(2.7.2) 


This  expression  is  similar  to  the  formula  giving  the  deflections 
of  columns. 

Equation  (2.7*2)  can  be  rewritten  in  the  form 

w 

o  ^crit  T  lo 
o 
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It  Is  clear  that  if  la  plotted  against  ,  the  result 

will  he  a  straight  line  of  slope  .  This  is  known  as  the  Southwell 

plot  and  it  provides  a  method  of  determining  T  ..  starting  from  exper- 

crio 

Imental  data  on  T  as  a  fvmction  of  w  (Fig. 32). 

0  o 

When  the  initial  and  the  buckled  shapes  are  different  (i.e., 
not  proportional  to  ^  )  the  relationship  between  w^  and  is 

different  and  can  be  derived  from  energy  considerations  (see  Ref.5)> 
one  obtains 


crlt 


1  + 


or 


w 

w  =  T  ..  -  K,w. 

o  crit  T.  1  io 
0 

where  is  a  coefficient  depending  only  on  the  initial  shape. 

Again  the  slope  of  the  straight  line  resulting  from  plotting  w^ 

versus  w  /t  determines  T  ... 
o'  o  crit 
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SECTION  III 


EXPERIMENTAL  METHODS  AND  RESULTS 


3.1  Introduction 

Instability  Is  not  always  characterized  by  a  sudden^  clearly  dis¬ 
tinguishable  deflection  from  the  original  configuration  but  rather,  In 
most  practical  cases,  by- an  Increase  In  the  rate  of  deflection.  This 
Is  due  to  Imperfections  In  either  the  Initial  configuration,  or  the 
application  of  the  load. 

The  determination  of  the  critical  time  of  instability  by  means  of 
direct  observation  is  possible  only  for  those  configurations  In  which 
the  ideal  conditions  are  closely  attained.  This  Is  not  the  case, 
unfortunately,  for  thin  plane  plates.  The  Initial  Imperfections  of 
manufacture,  the  effect  of  gravity,  the  perturbations  Introduced  by  the 
attachment  of  temperature  and  deflection  measuring  devices,  and  perhaps 
other  causes,  such  as  imperfections  in  the  supporting  rig,  make  initial 
flatness  unobtainable  in  experimental  work. 

The  Initial  configuration  is  not  flat  and  consequently  deflections 
appear  as  soon  as  the  heat  Is  applied  and  stresses  develop  In  the  plate. 
Nevertheless  the  critical  time  of  the  Ideal  case  can  be  determined  from 
experimental  data  If  the  methods  developed  in  Section  II  are  applied. 

The  critical  temperature  can  be  calculated  from  data  on  the  tenqper- 
ature  distribution  by  means  of  the  linear  theory  described  In  2.4.  This 
was  done  in  article  2.6  for  a  particular  case. 

Another  way  of  determining  the  critical  temperature  Is  by  means  of 
the  Soutlwell  plot  described  In  2.7.  This  method  requires  data  on  the 
relation  between  the  temperature  and  the  deflection  at  the  inner  edge 
of  the  plate  during  the  process  of  buckling. 

In  the  following  articles  a  summEury  Is  given  of  the  methods  applied 
or  developed  in  this  research  effort  for  the  gathering  of  data  on  teniper- 
ature  emd  deflections;  the  tests  performed  are  described  and  the  data 
obtained  are  used  to  determine  the  critical  tenqwrature  by  means  of  a 
Southwell  plot.  Finally  the  results  eure  conpared  with  the' predictions  of 
the  previously  obtained  formulas. 
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3.2  Measurement  of  the  Lateral  Deflections 


The  selection  of  a  displacement  measuring  device  for  a  particular 
experiment  should  depend  upon  a  number  of  considerations  relating  to 
the  specific  characteristics  of  that  experiment.  Perhaps  the  most 
important  consideration  is  the  order  of  magnitude  of  the  displacements 
to  be  measured. 

In  our  particular  case  we  are  Interested  in  corroborating  theoret¬ 
ical  results  obtained  from  the  linear  theory  on  the  deflection  of  thin 
plates.  Consequently  we  are  concerned  with  deflections  smaller  than 
the  thickness  of  the  plate. 

Since  the  plate  specimens  tested  in  our  experiments  had  a  thickness 
of  h  =  0.058  in.,  we  required  a  device  capable  of  measuring  with  rea¬ 
sonable  accuracy  displacements  in  the  range  from  0  to  50  thousandths  of 
an  inch. 

The  criterion  of  what  can  be  considered  a  reasonable  accuracy  is 
the  degree  of  precision  required  in  the  determination  of  the  critical 
temperature.  In  an  average  test  it  was  found  that  an  absolute  error 
in  the  temperature  reading  of  ±  2°F  and  an  absolute  error  in  the 
central  deflection  residing  of  0.001  in.  result  in  a  relative  error 
of  8  per  cent  in  the  value  of  the  critical  temperature  as  calculated  by 
the  Southwell  plot.  This  is  considered  acceptable. 

Thus  the  requirements  are  a  usable  rsmge  from  0  to  50  thousandths 
of  an  inch  with  an  absolute  error  smaller  than  one  thousandth  of  an  inch. 

These  are  not  hard  requirements  to  fulfill  and,  as  a  natter  of  fact, 
a  regular  dial-gauge  satisfies  them.  Nevertheless,  other  considerations 
rule  out  the  use  of  dial  gauges  in  favor  of  an  optical  system. 

The  advemtages  of  an  optical  system  are  veu'Ious;  two  of  them  are 
of  peurtlcular  loqportance: 

(a)  An  optical  system  does  not  Interfere  with  the  phenomena,  emd 

(b)  it  provides  a  direct  picture  of  the  deflection  pattern  all 
over  the  plate. 
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The  optical  system  of  deflection  measurement  selected  for  our 
experiments  Is  schematically  shown  In  Fig. 12.  A  parallel  beam  of  rays 
of  light  carrying  the  Image  of  a  grid  of  straight  lines  parallel  to  the 
plane  xz  (the  plane  xz  Is  determined  by  the  normal  to  the  plate  and 
the  direction  of  incident  light)  is  projected  on  the  surface  of  the 
specimen,  which  has  a  mirror  finish,  and  Is  reflected  by  the  specimen 
on  an  observation  milk  glass  G  parallel  to  the  plane  of  the  specimen. 
Thus  an  Image  is  formed  on  the  milk  glass. 

When  the  plate  Is  flat  the  Image  of  the  straight  line  abc  Is  a 
straight  line  a^b^c^  .  If  it  Is  assumed  that  the  plate  Is  fixed  at 
Its  outer  edge,  the  slopes  of  the  elements  a'  and  c'  will  not  change 
when  the  plate  deflects  and  consequently  their  Images  a^  and  c^  will 
not  move. 

Let  us  consider  now  the  element  b'  on  the  axis  Oy  .  The  normal 
to  this  element,  in  the  Initial  flat  state,  is  n^  In  the  plane  yz  . 

If  the  deflections  of  the  plate  are  axisymmetrlcal,  the  element  b' 
will  deflect  to  a  position  such  as  bg  and  its  normal  will  be  Ug  in 
the  plane  yz  .  Thus  the  angle  (p  between  n^  and  ng  is  equal  to 
the  change  in  radial  slope  of  the  element  b’  and  the  image  of  bg 
will  be  a  certain  point  bg  . 

If  the  deflection  w  of  the  element  b'  is  neglected  conpeured 
with  the  distance  d  between  the  planes  of  the  plate  imd  the  obser¬ 
vation  glass,  it  follows  from  elementary  optics  and  geometry  that 

tan  (p  Z!  (p  =  ^  (3.2.1) 

where  D  is  the  distance  b^^bg  ,  that  is  the  central  distortion  In 
the  direction  parallel  to  Oy  of  the  image  of  the  line  abc  . 

The  photographs  of  Figs.  13  and  14  show  how  the  IntEige  of  the  grid 
of  lines  appeeurs  for  the  plate  in  the  flat  state  and  in  the  deflected 
state.  The  Images  corresponding  to  an  intproved  grid  eure  given  in 
Fig. 15. 
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By  comparison  of  two  photographs  corresponding  to  the  Initial  and 
the  deflected  states  the  Image  displacements  D  of  a  number  of  points 
on  the  diameter.  Oy  can  he  measured  and  the  corresponding  changes  In 
slope  can  be  calculated.  Thus  the  function  ^(r)  =  dw/dr  can  be 
plotted  and  by  Integration  the  deflection  of  any  point  will  be  given 
by 


(3.2.2) 


3.2.1  The  Moire  Technique 

The  Image  distortion  D  can  be  magnified,  If  necesseury,  by 
means  of  the  so-called  Moire  fringe  technique.  This  technique  Is 
widely  known  and  has  been  the  subject  of  a  number  of  papers  some  of 
which  are  listed  under  the  references. 

Basically  the  method  consists  In  superlnqioslng  on  the 
Image  of  a  grid  of  lines  (usually  parallel  lines  either  marked  on  the 
surface  of  the  model  to  be  studied  or  reflected  by  Its  surface  on  an 
observation  screen,  as  In  our  case)  another  analogous  grid  of  lines, 
slightly  rotated  with  respect  to  the  former  one,  so  eis  to  produce  a 
so-called  mechanical  Interference  resulting  In  fringes  which  look  some¬ 
what  like  the  textile  material  known  as  Moire  (from  here  the  name  of 
the  phenomenon) . 

When  the  model  Is  unloaded  the  Image  of  the  projected  grid 
of  lines  appears  as  an  array  of  parallel  stralf^t  lines  (Flg.l6),  and 
superposition  on  this  Image  of  an  analogous  grid,  sllj^tly  rotated 
(Fig. 17)  results  In  a  Moire  pattern  with  stralf^t  fringes  (Flg.lS). 

When  the  model  Is  loaded  the  Image  of  the  projected  grid 
appears  distorted  (Flg.l9)  and  superi>osltlon  of  the  Interfering  grid 
results  In  a  Moire  pattern  of  ciurved  fringes  (Fig. 20). 
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If  the  relative  rotation  of  the  two  grids  Is  snail,  small 
distortions  of  the  projected  image  result  in  large  distortions  of  the 
Moire  pattern. 

In  the  experiments  conducted  in  this  research  effort  the 
distortions  D  of  the  projected  image  were  large  enough  to  be  measured 
directly  and  consequently  no  use  was  made  of  the  Ktoire  technique  in 
the  calcxxlatlon  of  the  deflections.  The  Interpretation  of  the  distor¬ 
tion  of  the  Moire  pattern  is  included  in  Appendix  A.  It  was  felt  that 
the  complications  Involved  in  the  calculation  of  the  distortions  D 
by  interpretation  of  the  Moire  pattern  were  not  Justified,  as  a  direct 
reading  of  D  from  an  enlarged  view  of  the  image  obtained  in  an 
Optical  Comparator  was  sirfficlently  accurate.  Nevertheless,  from  the 
qualitative  point  of  view,  the  Moire  technique  proved  to  be  very 
satisfactory  for  direct  observation  of  Imperfections  in  the  initial 
flatness  of  the  plates,  as  well  as  for  the  observation  of  the  rate  of 
deflection  during  the  buckling  phenomenon. 

Thus  it  was  decided  to  use  the  Moire  technique  in  half  of 
the  image,  for  qualitative  observation,  leaving  the  other  half  as  direct 
image,  for  quantitative  determination  of  the  deflections. 

The  resulting  comrposlte  image  is  shown  in  Figs. 21  and  22. 
Figure  21  corresponds  to  the  initial  state,  the  upper  half  being  a 
direct  image  as  reflected  by  the  plate  and  the  lower  part  a  fringe 
pattern  obtained  by  the  Moire  technique.  Fig. 22  is  the  image  in  the 
buckled  state. 

The  semicircles  appeeurlng  on  the  direct  image  seirve  to  check 
the  axisymmetry  of  the  deflected  mode. 

3.3  Description  of  the  Test  Rig 

The  test  rig  is  sketched  in  Fig. 23  and  photographed  in  Figs. 24  and 
25.  It  consists  of  an  optical  bench  on  which  a  300  watt  concentrated 
arc  leunp  provides  a  light  source  of  high  brlllleuice  and  small  size. 

The  light,  after  being  projected  through  the  target,  is  converted  by 
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means  of  a  suitable  set  of  lenses  Into  a  parallel  beeun  of  light  large 
enough  to  cover  the  whole  surface  of  the  plate  to  be  tested.  A  mirror 
reflects  the  beam  of  llgd^t  on  the  plate  which  Is  movinted  horizontally 
on  a  massive  support  (Flg.26).  The  highly  polished  surface  of  the 
plate  reflects  the  light  onto  an  observation  milk  glass  where  the  Image 
Is  formed.  The  observation  mirror  situated  above  the  milk  glass  Is 
only  for  convenience  of  observation. 

The  temperature  distribution  over  the  plate  Is  recorded  by  means 
of  chrome 1-alumel  thermocouples  attached  to  the  lower  surface  of  the 
specimen.  The  axlsymmetrlc  character  of  the  heating  Is  also  checked. 

The  heating  unit  (Flg.27)  was  designed  so  as  to  approach  as  closely 
as  possible  the  conditions  of  radial  heating.  The  heater  consists  of 
a  resistance  wire  element  laid  on  a  support  of  asbestos.  When  mounted 
on  the  plate  Its  two  halves  clamp  on  the  Inner  edge  of  the  plate, 
slightly  overlapping  it.  The  heater  weighs  only  two  ounces  and  moves 
with  the  plate  during  the  process  of  buckling. 

The  specimens  used  In  the  tests  were  annular  plates  with  the 
following  characteristics: 

Material:  steel 

Poisson  ratio:  v  =  0.25 

Coefficient  of  thermal  expeuision:  a  »  6.1xl0”^  °f”^ 

Inner  radius:  a  =  1.5  in. 

Outer  effective  radius:  b  =  5  in. 

Thickness:  h  =  0,058  In. 

The  front  surface  of  the  specimens  was  polished  to  a  mirror 
finish. 

3.1+  Buckling  Tests 

Three  tests  were  conducted  for  different  values  of  the  heat  flux. 
The  changes  In  the  beat  flux  were  obtained  by  means  of  a  variable 
rheostat  In  series  with  the  heating  imlt. 
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A  continuous  record  of  the  temperature  distribution  vas  obtained 
and  a  simultaneous  sequence  of  pictures  ws  talten  of  the  image  diirlng 
the  heating  process. 

The  sequence  of  images  obtained  during  test  no.l  is  shown  in  Fig. 28. 
The  photographs  correspond  to  0,  30,  6l,  90,  120,  150  and  l80  seconds 
respectively. 

By  comparing  these  photographs  with  a  pictxire  of  the  image  as 
reflected  by  a  calibration  flat  mirror,  the  displacements  D  (propor¬ 
tional  to  the  change  in  slope)  can  be  determined.  This  can  be  done 
with  great  accuracy  in  an  optical  conqmrator. 

In  this  manner  the  change  in  slope  distribution  over  the  radius 
can  be  plotted  for  the  recorded  Intervals  of  time.  This  has  been  done 
in  Fig. 29  for  test  no.l.  In  this  plot  B  is  the  slope  in  the  radial 
direction  and  K  is  a  constant  which  Includes  the  factor  l/2d  (see 
Eq. (3.2.1))  and  the  magnification  factor  of  the  optical  conqparator. 
Integration  of  these  curves  yields  the  deflections. 

As  mentioned  before,  in  this  particular  test  the  Images  were  com¬ 
pered  with  those  of  a  flat  calibration  mirror.  Consequently,  the  slopes 
rather  than  the  cheinges  in  slope,  were  obtained  for  every  point.  It 
cem  be  seen  that  for  t  =  0  the  central  deflection  appears  to  be 
Wj^^  =  3.76  thousandths  of  an  inch. 

The  verlatlon  of  the  total  central  deflection  w.  +  w  with  time 

10  o 

is  shown  in  Fig. 30.  Actually,  for  the  purpose  of  determining  the 
critical  temperature  by  means  of  the  Southwell  plot,  the  determination 
w^^  is  not  required.  In  tests  2  and  3,  the  images  of  the  plate  for 
the  successive  intervals  were  compared  with  the  image  of  the  initial 
state  (  t  =  0  ).  In  this  manner  the  changes  in  slope  eoe  directly 
determined  and  the  corresponding  Increments  in  the  deflection  can  be 
calculated. 

In  Fig. 31,  corresponding  to  test  no.l,  the  edge  tenqperature 

Increment  T  cuid  the  central  deflection  Increment  w  eure  plotted 
o  o 

against  time.  From  these  data  a  Southwell  plot  cem  be  made.  This  is 
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done  In  Fig. 32  where  one  can  see  that  the  Southwell  line  la  straight  In 
the  small  deflection  region.  It  determines  the  value 


(’'crjj  ■ 

Tests  2  and  3  were  conducted  In  the  sane  manner  (see  Figs. 33,  34, 
35,  36  and  37)  resulting  In  the  following  critical  tenqperatures ; 

(''crlt)g  ' 

('crlt),  ■ 

For  the  purpose  of  comparison,  these  results  were  plotted  In  Fig. 11 
which  summarizes  the  results  of  hoth  theory  and  experiment. 
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APPENDIX  A 


INTERPRETATION  OF  THE  DISTORTIONS  IN  THE  INTERFERENCE 
FRIN®S  OBTAINED  BY  THE  MOIR^  TECHNIQUE 

The  nature  of  the  Moire  interference  pattern  can  be  better  under¬ 
stood  by  observation  of  the  phenomenon  on  a  magnified  scale.  In  Fig.A-1 
it  can  be  seen  that  the  direction  and  location  of  the  black  interference 
fringes  correspond  to  those  of  lines  such  as  a  along  which  there  is 
total  coverage  of  light. 

If  the  interfering  grids  are  schematically  represented  by  thin 
lines  (Fig.A-2),  the  location  of  the  black  Interference  fringes  is 
determined  by  the  middle  points  of  the  segments  between  intersection 
points. 

Using  this  schematic  representation  Flg.A-3  has  been  constructed. 

Let  a^  and  a^  be  two  consecutive  lines  of  the  image  projected 
by  the  specimen  in  the  initial  state ^  and  let  b  be  one  line  of  the 
fixed  Interfering  grid. 

If  A  is  the  middle  point  of  the  segment  determined  by  the  inter¬ 
section  of  b  with  a^  and  a^  ,  there  will  be,  in  the  initial  state, 
a  certain  fringe  f  passing  through  the  point  A  . 

When  the  specimen  deflects,  the  lines  a^  emd  a^  become  a^ 
and  ag  and  the  point  A  moves  along  the  line  b  to  the  point  A', 
the  fringe  f  becoming  f  . 

The  displacement  D  of  the  point  P  can  then  be  obtained  from 
reading  the  distance  AA'  =  6 

D  =  5  sin  a  6a 

where  a  is  the  euigle  formed  by  the  grids. 

Nevertheless,  as  it  was  discussed  in  3.2,  the  displacement  D  is 
proportional  to  the  change  in  radial  slope  only  when  P  is  on  the 
diameter  Oy  euid  the  deflection  is  oxisymmetrlc. 
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The  displacement  D  of  points  P  that  are  not  on  the  axis  Oy 
determines  only  one  component  of  the  maximum  slope  of  the  element 
corresponding  to  point  P  . 

To  obtain  readings  of  D  corresponding  to  points  P  on  the  axis 
Oy  only  the  points  A*  lying  on  Oy  can  be  considered. 

In  our  case  It  was  found  that  the  niimber  of  fringes  that  Inter¬ 
sected  the  axis  Oy  was  not  sxifflclent  to  obtain  a  satisfactory  deter¬ 
mination  of  the  distribution  of  change  In  slope  along  Oy  . 

In  conclusion  It  was  decided  that  the  Moire  technique  was  satis¬ 
factory  only  for  qualitative  observation  while  the  method  of  direct 
Image  as  explained  In  3,2  was  more  reliable  and  useful  for  the  quanti¬ 
tative  estimation  of  the  changes  In  slope  of  the  plate. 
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Fig.l.  tteasperature  Distribution  from  Theory  Assuming 
no  Heat  Dissipation  to  the  Environment 
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Fig. 4.  Test  Hb.l  -  Teagierature  Distribution  over  the  Plate 
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Fig. 7,  Stodola-Vianello  Method.  Graphical  Integration 
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Fig. 8.  Edge  Temperature  Variation  with  Time  (Theory  -  no  Heat  Dissipation) 


2 

Fig.  9*  Critical  Edge  Tenperature  (T  /h  )  and  Critical  Time 

o  crit 

versus  Heat  Input  Flux  (Theory  -  no  Heat  Dissipation) 


Fig.l4.  Image  of  Loaded  Plate 
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Fig. 15.  (a)  Image  of  Unloaded  Plate 


Fig. 15.  (b)  Image  of  Loaded  Plate 
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Fig.18.  Moire  Fringes  of  Fig. 20.  Moire  Fringes  of 

Unloaded  Plate  Loaded  Plate 
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Fig. 21.  Image  of  Unloaded  Plate 


Fig. 22.  Image  of  Loaded  Plate 
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Flg.23.  Test  Rig 


■7 


Flg.2U.  Test  Rig 


Fig. 26.  Plate  Mounted  on  its  Support 


Fig. 27.  Plate,  Support  and  Heater 
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Fig. 28.  Test  No.l  -  Sequence  of  Images 
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Fig. 29.  Test  No.l  -  Change  In  Slope  Distribution  with  Time 
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Fig.  31.  Test  No.l  -  Central  Edge  Temperature  and  Central 

Ed^  Deflection  versus  Tine 
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Fig. 33.  Test  No. 2  -  Change  In  Slope  Distribution  with  Time 
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Fig. 34.  Test  Ho. 2  •  Centred.  Edge  Temperature  and  Central 

Edge  Deflection  versus  Time 
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Fig. 36.  Test  No. 3  -  Central  Edge  Temperature  and  Central 

Edge  Deflection  versus  Time 
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Fig. 37.  Test  No. 3  -  Southwe31  Plot 
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